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We use a series of ray-tracing experiments to determine the magnification distribution of high-redshift sources 
by gravitational lensing. We determine empirically the relation between magnification and redshift, for various 
cosmological models. We then use this relation to estimate the effect of lensing on the determination of the 
cosmological parameters from observations of high -2 supernovae. We found that, for supernovae at redshifts 
2 < 1.8, the effect of lensing is negligible compared to the intrinsic uncertainty in the measurements. Using 
mock data in the range 1.8 < 2 < 8, we show that the effect of lensing can become significant. Hence, if a 
population of very-high -2 supernovae was ever discovered, it would be crucial to fully understand the effect of 
lensing, before these SNe could be used to constrain cosmological models. 


1. INTRODUCTION 

High-redshift supernovae have become a major tool 
in modern cosmology. By measuring their appar¬ 
ent magnitude^ we can estimate their luminosity dis¬ 
tances di (see [lllallSj references therein). Since 
the relationship between dr and the redshift z de¬ 
pends on the cosmological parameters, observations 
of distant SNe can constrain the cosmological model. 
Prior to the announcement of the WMAP results Q, 
observations of high-z SNe provided the most com¬ 
pelling evidence of the existence of a nonzero cosmo¬ 
logical constant. 

The luminosity distances d^ are determined by 
combining the observed fluxes F with estimates of the 
SNe luminosities L. Uncertainties in are caused by 
uncertainties in L, because SNe are not perfect stan¬ 
dard candles. The flux F is much easier to measure, 
but for distant sources the value of F might be al¬ 
tered by gravitational lensing cause by the interven¬ 
ing distribution of matter. For instance, a positive 
magnification would result in a increase in F, and an 
underestimation of dr- 


2. THE ALGORITHM 

We have developed a multiple lens-plane algorithm 
to stu^ light propagation in inhomogeneous universes 
[M Q algorithm, the space between the 

observer and the sources is divided into a series of 
cubic boxes of comoving size 128 Mpc, and the matter 
content of each box is projected onto a plane normal 
to the line of sight. The trajectories of light rays are 
then computed by adding successively the deflections 
caused by each plane. 

To use this algorithm, we need to provide a descrip¬ 
tion of the matter distribution along the line of sight. 
Matter is divided into two components: background 
matter and galaxies. We use a P^M algorithm to sim¬ 
ulate the distribution of background matter. The sim¬ 


ulations used 64^ equal-mass particles and a 128^ PM 
grid, inside a comoving volume of size 128 Mpc. The 
matter distribution in the different cubes along the 
line of sight then corresponds to the state of the sim¬ 
ulation at different redshifts.^ We then use a Monte 
Carlo method for locating galaxies into the compu¬ 
tational volume 0,1^ . Galaxies are located accord¬ 
ing to the underlying distribution of background mat¬ 
ter. Morphological types are ascribed according to the 
morphology-density relation 0. Galaxies are mod¬ 
eled as nonsingular isothermal spheres, with rotation 
velocities and core radii that vary with luminosity and 
morphological types. In Figure Q we use a series of 
zooms to illustrates the dynamical range of the den¬ 
sity distribution generated by this method. 


3. THE RAY-TRACING EXPERIMENTS 

We consider 3 cosmological models: (1) a flat, cos¬ 
mological constant model with Hq = 0.27, Aq = 0.73, 
and Hq = 71km/s/Mpc. This model is in agree¬ 
ment with the results of the WMAP satellite 01 ■ (2) 
a low-density model with Hq = 0.3, Aq = 0, and 
Hq = 75km/s/Mpc. (3) a matter-dominated model 
with Ho = 1) Ao = 0, and Hq — 75km/s/Mpc. For 
each model, we consider sources at 5 different red- 
shifts: Zs = 1, 2, 3, 4, and 5. For each combination 
mo del-redshift, we performed 10-20 ray tracing ex¬ 
periments. Each experiment consists of propagating a 
square beam of 101 x 101 = 10, 201 rays back in time 
from the present to redshift Zg, through the matter 
distribution. The rays in the beam were widely sepa¬ 
rated, by 6 arc minutes, and therefore sampled differ¬ 
ent regions of space. We computed the magnification 
matrix A along each ray. The magnification g is then 


^In practice, we combine cubes from different simulations in 
order to avoid periodicities along the line of sight. 
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given by 


1 

det A 


( 1 ) 


Figure [3 shows the distribution of magnifications for 
the A-model. The distribution peaks at /i < 1, and 
is strongly skewed. The width of the distribution in¬ 
creases with the source redshift. The distributions for 
the other two models are qualitatively similar. 


4. THE EFFECT OF LENSING ON 
STATISTICS OF HIGH-Z SUPERNOVAE 


Estimating the effect of lensing on the statistics of 
high-z supernovae is a complex problem, and we defer 
such analysis to future work. Here we use a simple ap¬ 
proach. First, for each model and each source redshift, 
we compute the standard deviation cr^ of the magni¬ 
fication distribution The values are shown in 

Figure 13 We use an empirical fit of the form 


(T 


r- ~ 


az 

l + bz ■ 


( 2 ) 


This enables us to estimate the values of cr^ at any 
redshift. 

We estimate the effect of lensing as follows: the 
distances of high -2 supernovae are reported in the lit¬ 
erature as: 


log(dLi?o) = a±(5a, (3) 


where di, is the luminosity distance, idg is the Hubble 
constant, a is the measurement, and 5a is the intrinsic 
uncertainty (i.e. not caused by lensing). The distance 
dL is related to the luminosity L and flux F by 


F =- 

Airdl 

After eliminating d^, we get 


( 4 ) 

( 5 ) 


The effect of lensing will be to modify the flux F. To 
account for it, we replace F by FFAF in equation ©, 
and expand to first order in da and AF. After some 
algebra, we get 

log(dLHo) =a±Sa± ■ (6) 

The last term represents the effect of lensing. We then 
make the approximation AF/F 'zz and get 

log(dLido) = a ± do ± , (7) 

where 5^{z) = (T^( 2 )/ 21 nl 0 is computed using the 
empirical relations plotted in Figure |3 We use the 


values of a and do reported by Tonry et al. M (their 
Table 8). In Figure 0 we plot the ratio d^/do versus 2 
(left of the dashed line). This quantity increases with 
redshift, but never gets higher than 0.5 for the Tonry 
et al. sample. Furthermore, we shall assume that 
do and d^ are statistically independent, and combine 
them in quadrature, using 

5 = {Sl + 5l)^/\ ( 8 ) 

where d is the total error. The contribution of lensing 
to this error is then of order 25% at most. 

The top panel of Figured shows a Hubble diagram 
[deviation A(m — M) relative to an empty universe, 
versus redshift], obtained by averaging the data in red¬ 
shift bins, using 

Wi = l/5f, (9) 

[A(m - M)]^. = EiWiA{m - M)/T,iWi, (10) 

6 , = (l/E.zc,)'/", (11) 

where the sums are over all data points i in bin j. 

The three curves, from top to bottom, show the ex¬ 
act results for the ACDM, low-density, and matter- 
dominated models, respectively. The results support 
the flat ACDM model and exclude the other models 
considered. 

The other panels of Figure |3 show the effect of 
lensing (the three models have to be plotted sepa¬ 
rately, because the correction due to lensing is model- 
dependent). This effect is totally negligible. The 
largest correction to the error bars is about 10% for 
the highest redshift bin, for the ACDM model. 

Clearly, the potential error introduced by lensing is 
negligible in comparison to the intrinsic error in the 
measurement, at least for SNe at redshifts 2 < 1.8. 
The next step is to estimate the effect on a yet- 
undiscovered population of very-high -2 SNe. We gen¬ 
erated a mock catalog of 43 SNe in the range 1.8 < 2 < 
8. We assume that the values of a are consistent with 
a ACDM model, and the values of 5a in that range are 
comparable to the ones in the range 1.5 < 2 < 1.8. 

The ratios d/j^/da are plotted in Figure ^ The ef¬ 
fect of lensing rapidly becomes important at redshift 
2 > 2. Figure El shows the Hubble diagrams of Fig¬ 
ure El which have been extended to higher redshifts to 
include the mock data. The error bars get significantly 
bigger when lensing is included. Furthermore, at red¬ 
shift 2 « 3, it becomes very difficult to distinguish 
the open, low density model from the cosmological 
constant model. 


5. SUMMARY AND CONCLUSION 

We have performed a series of ray-tracing experi¬ 
ments using a multiple lens-plane algorithm. We have 
estimated the standard deviation of the magnification 
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distribution with source redshift, for three different 
cosmological models. Using this relation, we have es¬ 
timated the effect of leasing on the statistics of high- 
redshift supernovae. The errors introduced by lensing 
are unimportant for SNe with redshift z < 1.8. How¬ 
ever, the effect of lensing on a hypothetical population 
of SNe at redshifts z > 1.8 could be very important, 
and must be understood before these SNe could be 
used to constrain cosmological models. 
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Figure 2: Distribution of magnifications for the fio = 0.27, Ao = 0.73 model. The various curves correspond to different 
source redshifts Zs, as labelled. 



z 


Figure 3: Standard deviation versus redshift, for all three models considered. The solid lines show empirical fits of 
the form cr^ = az/{l + bz). 
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Figure 5: Hubble diagram showing the magnitude deviation A(m — M) relative to an empty universe, for the three 
models considered. In the top panel, the three curves, from top to bottom, show the analytical result for the 
cosmological models (f2o,Ao) = (0.27,0.73), (0.3,0.0), and (1.0,0.0), respectively. The last three panels reproduce the data 
of the top panel, but have been corrected to account for lensing. Since this correction is model-dependent, the three 
models are plotted on separate panels. Error bars show 90% confidence level. 
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Figure 6: Same as Figure El but the panels have been extended to include the mock data (four rightmost data points). 


1415 




























